Submitted to ApJL on June 2, 2003 

Preprint typeset using I4T]^X style emulateapj v. 3/16/03 



m 
O 
o 

(N 
C 

(N 



NEUTRON-RICH HYDROMAGNETIC OUTFLOWS IN GAMMA-RAY BURST SOURCES 

Nektarios Vlahakis, Fang Peng, and Arieh Konigl 

Dept. of Astronomy & Astrophysics and Enrico Fermi Inst., Univ. of Chicago, 5640 S. Ellis Ave., Chicago, IL 60637 
vlahakis@jets.uchicago.edu, fpeng@oddjob.uchicago.edu, arieh@jets.uchicago.edu 
Submitted to ApJL on June 2, 2003 

ABSTRACT 

We demonstrate that "hot" MHD outflows from neutron-rich black-hole debris disks can significantly alleviate 
the baryon-loading problem in gamma-ray burst (GRB) sources. We argue that the neutron-to-proton ratio in 
disk-fed outflows might be as high as ~ 30 and show, with the help of an exact semianalytic relativistic-MHD 
solution, that the neutrons can decouple at a Lorentz factor 7d ~ 15 even as the protons continue to accelerate 
to 7oo ~ 200 and end up acquiring ~ 30% of the injected energy. We clarify the crucial role that the magnetic 
field plays in this process and prove that purely hydrodynamic outflows must have 7d > a few xlO 2 . The 
motion of the decoupled neutrons is not collinear with that of the decoupled protons, so, in contrast to previous 
suggestions based on purely hydrodynamic models, the two particle groups do not collide after decoupling. If 
the decoupled neutrons move at an angle > 1 w 3.8°(15/7d) to the line of sight to the GRB source, most of 
their emission after they decay into protons will remain unobservable. 

Subject headings: accretion disks — gamma rays: bursts — MHD — nucleosynthesis — relativity 
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1. INTRODUCTION 

Gamma-ray burst (GRBs) are inferred to arise in rela- 
tivistic outflows of te rminal Lorentz factor 7^ > 10 2 (e.g., 
iLithwick & Sarill2001l) . In the case of long-duration (At > 

2 s) bursts, there is growing evidence that the outflows are 
highly collimated and involve a k inetic energy ~ 10 51 ergs 
(e.g. JPanaitescu & Kumarll2002l) . The outflows likely origi- 
nate in a newly formed neutron star (NS) or stellar-mass black 
hole (BH) and are powered either by the NS or BH spin or by 
the gravitational potential energy of a remnant accretion disk. 
The jets can be driven either thermally or magnetically. Mag- 
netic fields offer a natural means of tapping the rotational en- 
ergy of the source and of guiding, accelerating, and collimat- 
ing the flow. In a recent series of papers, Vlahakis & Konigl 
[2001, 2003a (VK03a), 2003b (VK03b)] summarized the ev- 
idence in favor of magnetic driving as well as previous work 
on this topic, and presented exact semianalytic solutions of the 
"hot" relativistic MHD equations demonstrating that Poynting 
flux-dominated jets can transform > 50% of their magnetic 
energy into kinetic energy of 7^ ~ 10 2 - 10 3 baryons. 

The baryonic mass involved in long-duration 
GRB outflows is estimated to be M p = Ek/^ooC 1 ~ 

3 x 10" 6 (£jfc/10 51 ergs)(7 oo /200) _1 M Q . This estimate 
poses an apparent difficulty for realistic source models. For 
example, a BH debris disk would need to be a factor ~ 10 4 
more massive even if 10% of its gravitational potential energy 
could be converted into outflow kinetic energy. A question 
thus arises as to how the baryon loading of the jets remains 
so low even as the energy deposition in the flow is highly 
efficient. A promising proposal to al l eviate this problem was 
made by Fuller, Pruet, & Abazaiian (2000), who suggested 
that, if the source of the outflow is neutron-rich, then the 
neutrons (which are only very weakly affected by the elec- 
tromagnetic field and are accelerated primarily by collisional 
drag with ions) could in principle decouple from the flow 
before the protons attain their terminal Lorentz factor. If 
the mass source for the flow is a NS then it will clearly be 
neutron-rich, and this is evidently the case also for outflows 
fed by a BH debris disk (see §[2j. However, it turns out that 
the decoupling Lorentz factor 7d in a thermally driven, purely 



hydrodynamic (HD) out flow is of the order of the inferred 
value of 7no (e.g.. iDerishev. Kocharovskv. & Kocharovskvl 
119991 iBeloborodov 2003b; see which h as so far limited 
the practical implications of the Fuller et al. (2000) proposal. 
In this Letter we demonstrate that a hydromagnetic neutron- 
rich outflow can go a long way toward addressing the baryon 
"contamination" problem. In §|2]we consider the composition 
of GRB disks and construct an illustrative disk model in the 
context of the supranova scenario for GRBs. A semianalytic 
MHD solution representing an initially neutron-rich outflow 
that undergoes decoupling at 7d <C 7oo is presented in § [3] 
We summarize and discuss the implications of our results in 

§EI 

2. SOURCE COMPOSITION 

iPruet. Wooslev. & Hoffmanl J20031) reviewed the various 
models that attribute GRB outflows to disks around newly 
formed stellar-mass BHs. They considered disk accretion 
rates ffidisk = [Musk/C^o s -1 )] = 0.01 - 1 and viscosity pa- 
rameters a = 0.01 —0.1, and calculated typical midplane elec- 
tron fractions Y e = ^2j(ZjXj/Aj) (where Zj, Xj, and Aj are the 
proton number, mass fraction, and atomic mass number, re- 
spectively) near the event horizon of a 3M Q BH in the range 
0.045 - 0.527 (with Y e decreasing with decreasing a for a 
given flidisk and increasing with decreasing mdisk for a given 
a). It was argued that material reaching the base of the out- 
flow at the top of the disk would remain neutron-rich if the 
midpl ane value of F P wer e low. Similar results were presented 
bv lBeloborodovl ( l2003bl) . 

As a complement to this calculation, we now consider de- 
bris disks of the type likely to arise in the supranova sce- 
nario for GRBs fVietri & Stella 111998b . In this picture, a 
supernova explosion leads to the formation of a rapidly ro- 
tating, supramassive NS (SMNS), which, after losing ro- 
tational support through a combination of electromagnetic 
and gravitational radiation on a timescale of weeks to years, 
eventually collapses to a black hole (in the process trigger- 
ing the GRB). In the original proposal, the debris disk was 
envisioned to represent the outer layers of an SMNS that 
were left behind when its central regions collapsed. How- 
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ever, IShibata. Baumgarte. & Shapirol (|2000) argued, on the 
basis of general-relativistic, 3D hydrodynamic simulations 
of uniformly rotating stars with a comparatively stiff equa- 
tion of state (EOS), that the collapse likely involves the en- 
tire SMNS. An alternative possibility of forming a disk in 
stars with a stiff EOS is for matter to be shed centrifu- 
gally in the equatorial plane as the SMNS evolves (through 
the loss of angular momentum) toward the gravitational- 
instability point. We illustrate this possibility by considering 
the evolution of such an SMNS along the mass-shedding limit 
from the maximum-mass point to the maximum-angular- 
velocit y (f2) point, where it becomes unstable to collapse 
(see [Cook. Shapiro. & Teukolskvl 119941). Using the EOS- 
M equilibrium model results of Cook et al. ( 1994), we infer 
that the newly formed BH would have a (gravitational) mass 
M BH = 2.04M Q and (assuming that none of the shed mass be- 
came unbound) would be surrounded by a debris disk of mass 

Coo k et all 111994) showed that, prior to the formation of the 
BH, there always exists a region in the equatorial plane out- 
side the SMNS within which particles on circular orbits are 
unstable to radial perturbations. We thus expect a bona fide 
accretion disk to form only after the SMNS collapses, and we 
approximate its initial extent as stretching from the innermost 
stable circular orbit of the final BH (~ 4GM Bil /c 2 « 12 km) 
to the equatorial radius (/? e w 22 km) of the Maximum-M con- 
figuration (i.e., AC3i ~ 10 km, where G3 is the cylindrical ra- 
dius and the subscript i denotes the base of the flow). Most 
of the mass would be initially concentrated near the radius of 
the maximum-0 SMNS (~ 16 km, henceforth taken to be the 
typical value of G5j), so the disk binding energy can be esti- 
mated as E b w GM B HM disk /G3i w2 x 10 52 ergs. If we identify 
a characteristic burst duration of 10 s with the accretion time 
of the disk, we obtain a fiducial accretion r ate of = 0.006. 
This rate is below the range considered by P ruet etalJ (2003) 
even as the characteristic midplane density (> 10 13 g cm -3 ) is 
much higher than in any of their models. Nevertheless, this 
disk and the inner regions of most of the P ruet et alJ (120031) 
models share the property that their cooling is dominated by 
neutrino emission under marginally optically thin conditions. 
If the outflows from these disks are strongly magnetized, then 
much of the internal heating could be due to Ohmic dissi- 
pation of wound-up magnetic field even as the transport of 
angular momentum might be dominated by large-scale mag- 
netic stresses (see VK03b). If, however, we represent both of 
these processes in terms of an a viscosity, then a rough es- 
timate of the parameter a for the SMNS-spawned disk can 
be obtained by equating the accretion time to the viscous 
angular-mo mentum transport time in a neutrino-cooled disk 
model (e.g., iPopham. Wooslev & Frverll 19991) . which yields 
a=1.2xl0" 3 . 

We can estimate the density of the innermost mass 
shell of the maximum-M SMNS that ends up in the disk 
by setting M Hlsk = \4wlgp ,( p .)/Gmi - (2GM/R e c 2 )] 
( Lorentz. Ravenhall. & Pethick 1993), which yields 
Pa = 1.5 x 10 32 dynes cm " 2 . Based on the results of 
Douchin & Haensel (2001), the corresponding density is 
p s pa 6 x 10 13 g cm -3 , and the shell lies near the bottom of 
the inner crust and is made up of free neutrons with X a ~ 0.8 

1 We neglect any additional mass from the fallback of the original super- 
nova explosion that might surround the BH, as this material would not be as 
highly magnetized as the SMNS and would therefore be unlikely to partake 
in the GRB outflow. 



and of neutron-rich composite nuclei with Z ~ 40 and 
Z/A ~ 1 /6. The disk mass thus comes mainly from the NS 
inner crust; the contribution of the outer crust (~ 10~ 4 M Q ) 
is negligible. After being shed from the NS, the composite 
nuclei are dissociated and eventually end up releasing their 
excess neutrons (with the remainder consisting mostly of 
4 He). The exact composition of the disk-ejected material is 
determined by nuclear reactions within the disk (but below 
the freezeout point, where the nuclear reaction timescale 
comes to exceed the vertical travel time of the ejected gas) 
as well as by deuteri um an d helium production within th e 
flow (e.g., iLemoind 120021 iPruet. Guiles. & Fuller! [2002). 
However, collision with neutrons just prior to decoupling 
are expecte d to break a part both a p articles and deuterons 
(e.g., lPruet etal.1 120021 |Beloborodov 2003b). In the limit 
that only free nucleons remain at the time of decoupling, the 
value of Y e for the NS-shed material would be > 1/30 and 
the neutron (n) to proton (p) ratio would be (1— Y e )/Y e < 29. 
Although the n/p ratio at the base of the outflow would be 
lower than this value, we expect it to be high enough to 
eventually give rise to n-p drift speeds ~ c that would lead to 
decoupling and result in the collisional breakup of the bulk of 
the remaining composite nuclei. 2 In the illustrative outflow 
solution presented in § [5] we use a constant pre-decoupling 
n/p ratio of 30. 

3. FLOW ACCELERATION AND NEUTRON DECOUPLING 

The equations of motion for the neutron component of the 
outflow are given by 



(1) 



where pon is the rest-mass density of neutrons, t/ n K = 
(7„c , 7„ V D ) their four- velocity, and 

^ K = ^PopPOn(C/ K -C/ i r) (2) 

the neutron-proton collisional-drag force density. 

The rest of the fluid - consisting of protons (with rest-mass 
density pop), radiation (with pressure Pg), and e + e" pairs (with 
pressure Pm) - is moving with four-velocity U K = (7c, 7V). 
The equations of motion of the charged fluid component are 



POp U k U v + (Pr+Pm) [ 4^- + if 



: -^EM J 7 >(3) 



incorporating the p-n drag force density -T K and the elec- 
tromagnetic force density F£ M = [J ■ E/c , (J°E+ J x B)/c]. 
Here E , B is the electromagnetic field and J v = (j° , J) the as- 
sociated four-current (with J°/c representing the charge den- 
sity). As in VK03a, we assume that the photons possess a 
blackbody distribution (under optically thick conditions) and 
the pairs a Maxwellian distribution with zero chemical poten- 
tial. In the large-temperature limit (9 = k%T /m e c 2 ^> 1) the 
latter corresponds to a 4/3 polytrope, fjyi oc 9 4 oc Pr; in this 
case Pm/Pr = I8O/71- 4 « 1.85 is a constant, so the matter and 
radiation can be treated as a single fluid. 

2 In the "viscous heating and optically thin neutrino cooling" approxima- 
tion to our debris-disk model, the disk is vertically isothermal and has a tem- 
perature 7di sk = 1.8 X 10 10 K and a midplane density of 1.9 X 10 13 g cnr 3 
at tDj. For the estimated value of M p (see §0 and with X n /Xp = 30, the 
density at the bottom of the disk surface layer that ends up in the flow is 
Pdisk - 1-2 x 10" g cm -3 (corresponding to ~ 0.1% of Mji^ being ejected). 
These values of and are comparable to those used in the models of 
IPruet etalj|2003l) . 
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The system Q-Q, together with Maxwell's laws and the 
proton and neutron mass conservation relations can in princi- 
ple be solved to yield the velocities of the charged {U K ) and 
neutral (t/ n K ) components as well as the electromagnetic field 
(E , B) and the thermodynamic quantities (pop , Pon , ©)• How- 
ever, a simpler method that is sufficient for our purposes is 
described below. 

During the initial (pre-decoupling) outflow phase, the neu- 
trons and protons are well coupled and their velocities are 
comparable: £/ n K » U K . By eliminating the drag force between 
equations and (0 and setting £/* « U K , we end up with the 
same system of equations as in the pure proton/lepton case — 
we only need to replace p by p 0n + p 0p = p 0p ( 1 + pon/po P )- By 
applying the same r self-similar model as in VK03a, we can 
study the neutron-rich outflow and evaluate all the physical 
quantities, including U K . Next, equation Q with U* ~ U K 
(or eq. 0) yields the force T K . Finally, equation (0 gives 
the drift velocity, which can be written (with the help of the 
continuity equation) as 

U£ - U K « -2m p U u U«/<j {) cp 0p . (4) 
So long as the drift velocity is much smaller than the proton 
velocity, the neutrons remain well coupled to the protons. In 
the general (nonradial) case there are two possible modes of 
decoupling, which can be dubbed "along the poloidal field" 
and "in the transfield direction." The former happens when 
the relative velocity between protons and neutrons along the 
poloidal magnetic field, Vji — V^ii, grows to ~ Vji w c, and the 
latter when V n ± becomes ~ AC3cos^/7Td yn , where ■& is the 
jet opening half-angle, AGS cos ?9 is its width, and Td yn ~ z/"fc 
is the dynamical timescale. (The subscripts || , _L denote com- 
ponents of a vector along the poloidal magnetic field and in 
the transfield direction, respectively.) 
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4Pm)/poc 2 (dashed), and Poynting-to-mass flux ratio —(c/4tT)EB l j > /'ypQC 2 V^ 
(dot-dashed) as functions of height (z) along the innermost fieldline. Here 
Po - POn + Pop below the decoupling point and po = p 0p above it. The baryon 
mass-normalized quantities are discontinuous at the decoupling transition, 
where they increase by a factor 1 + po n / Pop- 

We model the post-decoupling evolution of the charged 
fluid component using the same r self-similar model used 
in the pre-decoupling phase (but replacing po by po p ). Fig- 
ures [0 and I2 show a solution representing an outflow from a 
disk with an inner radius r m = 1.2 x 10 6 cm and outer radius 
r out = 2r m . The initial conditions for the innermost fieldline 
in = r m , 8i = 80°) are: pon/Pop = 30, pon+Pop = 10 5 g cm" 3 , 8 = 
2.2, B p = 10 12 G, = -10 14 G, ■& = 55°, V p = 0.6916c, and 
= 0.35c. 3 These conditions correspond to a rough equipar- 
tition between comoving magnetic and thermal pressures — 

3 The value of the model parameter F is 1 .05 before decoupling and 0. 1 



(B 2 — £ ,2 )/87r(PR + f > M) = 1.0, a total outflowing baryon mass 
M b = 2 JfrypoV -dS At = 9.22 x 10- 5 (Af/10s)M Q , and total 
energy E\ = 3.62 x 10 51 (Af/10 s) ergs. The various compo- 
nents of the initial total energy are [electromagnetic, enthalpy 
(including rest energy), kinetic] = [5.615, 15.631, 0.583] x 
MbC 2 . During the pre-decoupling phase, part of the enthalpy is 
converted into baryon kinetic energy (resulting in an increase 
in 7) and another is converted into electromagnetic energy 
(resulting in an increase in the Poynting-to-mass flux ratio; 
see Fig. [Q. Above z ~ 10 9 cm the fluid is practically cold 
(Pr + Pm <S Poc 2 £ ~ 1) and the magnetic field is starting 
to accelerate the matter: in this regime the increasing 7 corre- 
sponds to a decreasing Poynting-to-matter flux ratio. 




z(cm) 



FIG. 2. — Proton-neutron drift velocity along the poloidal magnetic field 
(solid), in the transfield direction (dot-dashed), and in the comoving proton 
frame (dashed) for the outflow described in Fig.fTI 

Figure [2] shows the various components of the drift ve- 
locity. At z = Zd ~ 7 x 10 12 cm, Vji — V n || ~ c an d me neu " 
trons decouple from the protons. In this particular solution 
the "along the flow" decoupling happens first, although at 
the time of decoupling the neutrons also have a nonnegligi- 
ble outward-directed transfield velocity component (V n j^ ~ 
0.1c <IC V n ||). The latter component arises from the ongoing 
magnetic collimation of the jet: by the time the flow reaches 
Zd, 1? is already decreased to ~ 7°. 4 At the decoupling point 
7 = 7d ~ 15.3, and the energy of the charged fluid compo- 
nent is £ C)d = 1.16 x 10 51 (Ar/10 s) ergs = 32% £ [with the 
remaining 68% residing in the neutron kinetic and rest ener- 
gies, 7dMbC 2 /(l + pop/Pon)]- In units of the proton rest energy 
M p c 2 = MbC 2 / ( 1 + /Oon/po P X the various components of this en- 
ergy are [electromagnetic, enthalpy (including rest energy), 
kinetic] = [201.79, 1.22, 14.31] xM p c 2 . 

Above Zd we continue the integration of equation (0 with 
JF K = 0. Figure [l] shows that the acceleration in this regime 
is magnetic: the electromagnetic energy decreases and the 
Lorentz force accelerates the protons to a terminal Lorentz 
factor 700 = 200. The final proton kinetic energy is Ek = 
JooM p c 2 £ c .d, indicating that the acceleration efficiency in 
the post-decoupling region is w 100%. 

The above solution demonstrates that, in contrast to the HD 
case considered by previous authors, neutron decoupling in an 
MHD outflow can occur with 7d <C 700 ■ To understand why 

after decoupling. We use the form of the n-p collisional cross section given 
in lDerishev et all! 1993) . with cr = 3 X 10~ 26 cm 2 . 

4 The value of 1? for the proton component changes little after decoupling; 
to obtain the exact value of the angle between the two particle species at the 
time of decoupling, the system of equations HI l3l must be solved without 
the approximation (/" as U K . 
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magnetic acceleration is crucial to this outcome, consider ap- 
proximating the local scaling of the Lorentz factor with cylin- 
drical radius by 7 oc BES" 5 ' and the fieldline (or streamline) shape 
by r ~ z oc 03 1+<52 , where Si , S2 are positive constants. Under 
the approximation of exact mass conservation for the proton 
component, -fpo p oc 03~ 2 . At the decoupling point t co ii « Td yn , 
where r co n = m v / aocpo v . 5 Using the above scalings, we obtain 

7d = 7 1 (P0paWmp7 1 ) 5l/(Ml+1 " 52) and r * = n(7d/7.) (1+52)/51 ■ In 
the HD case, the energy equation £7 oc 07 = const and the 

polytropic relation 4 oc p^ 3 always imply a linear increase 
of 7 with 03. Thus, <5i = 1 in this case, and the minimum value 
of 7d, 7dHDmin = 7iG°0p,incrrj/»Jp7i) 1/3 , is attained in a coni- 
cal flow (62 = 0). For typical values of the source size (at 
least 10 6 cm) and initial proton density 6 (at least a few times 
10 2 g cm -3 ), 7dHDmin is a few times 10 2 . On the other hand, 
in MHD flows where a significant part of the enthalpy can be 
initially transferred into Poynting energy, the pre-decoupling 
acceleration is slower, corresponding to Si < 1. (In the so- 
lution presented here, Si ss 0.13 and 62 ~ 0.15.) This makes 
it possible for 7d to be <C 7dHDmin- The energy deposited in 
the Poynting flux is subsequently returned to the matter as ki- 
netic energy, enhancing the fraction of the total energy used 
to accelerate the proton component. 

4. SUMMARY AND IMPLICATIONS 

We argued that GRB outflows emanating from accreting 
debris disks around stellar-mass black holes could possess 
a high (< 30) pre-decoupling neutron-to-proton mass ratio, 
and we demonstrated semianalytically, by means of a "hot" 
relativistic-MHD solution, that n-p decoupling can occur at a 
Lorentz factor (7d ~ 15 in the presented example) that is sig- 
nificantly lower than the terminal 7 of the protons (700 = 200 
in our solution). We explained why, in contrast, 7d cannot be 
smaller than a few times 10 2 if the acceleration is purely hy- 
drodynamic. In the outflow solution that we constructed, the 
protons, while constituting only ~ 3% of the ejected mass, 
end up with a kinetic energy ~ 10 51 ergs that represents ~ 1/3 
of the injected energy (with the remainder going to neutron 
kinetic energy). 

This could significantly alleviate the baryon-loading prob- 
lem in GRB source models. For example, in the illustrative 
SMNS-shed disk model that we considered, the outflow rate 
represents ~ 0.1% of the mass inflow rate. If a fraction of 

5 A more exact condition for the decoupling comes from eq. |4j, which im- 
plies V|| -V„|| ~ c «=> 2(m f ,/(T p 0p )d'Y/dz ~ 1, or r dyn ~ (2rf In In z)r coll . 

6 As a rough estimate, pop.i ~ M p /4nrfcAt, with the ejected proton mass 



this order of the disk mass can be ejected (with X n /X p that be- 
comes ~ 30 just before decoupling) in a "hot" magnetized jet, 
then it may not be necessary to consider alternative sources 
of energy or mass for GRB outflows, for which the converse 
problem (how to avoid havin g too few b aryons) is often en- 
countered (e.g., Levinson & Eichler 2003). 

The decoupled neutrons decay into protons at a distance 
r p = l&cTp w4x 10 (7d/15) cm, where 773 k, 900 s is the 
comoving decay time. Nonthermal shock-induced emission 
from the resulting charged-particle shell may be hampered by 
the lack of a strong magnetic field unless signifi cant fie ld am- 
plification occurs in the shocks themselves (e.g., Beloborodov 
I2003al) or if the neutrons propagate within a highly magne- 
tized medium (as might happen in the supranova scenario; 
Koniri & Granot 2002). Internal collisions arising from the 
same velocity perturbations that are invoked to account for 
the observed 7-ray emission in the internal-shock scenario 
would typically occur on a scale <C rp and thus (in con- 
trast to the situation in the proton shell) would not give 
rise to radiative shocks. Emission might, however, be ex- 
pected from the reverse shock driven into the decaying neu- 
tron shell as it starts to decelerate. Most of this shell's ki- 
netic energy would, however, be radiated from the external 
shock that is driven into the ambient gas. The bulk of the 
radiation would be emitted at the distance fdecn where the 
mass swept up from the ambient medium comes to exceed 
~ 1/74 of the shell mass. This would occur on a timescale 
~ ^dec,n/27|c, which, for a uniform environment, is a factor 

- 10 3 (X n /30X p ) 1 / 3 (7 oo /137d) 7/3 larger than the correspond- 
ing timescale for the original proton shell. Most of this emis- 
sion would, however, remain undetectable if the angle be- 
tween the neutron-shell velocity and the line of sight exceeded 
~ l/7d ~ 3.8°(15/7d). It is in any case clear that the trans- 
verse n-p relative motion at the time of decoupling precludes 
any su bsequent interactio ns between the two shells (see, e.g., 
Pruet & Dalai 2002 and Beloborodov 2003a). A detailed con- 
sideration of the possible observational signatures of the neu- 
tron shell remains an interesting problem for future research. 

We thank Andrei Beloborodov, Gregory Cook, Jason Pruet, 
and Stuart Shapiro for helpful discussions and correspon- 
dence. This work was supported in part by NASA grant 
NAG5-12635. 

inferred from M p = Ei i j r y aa c 1 using ~ 10 51 ergs and 700 = a few X 10 2 . 
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